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1  Phase  I  Technical  objectives  and  accomplishments 

In  Phase  I  we  developed  an  extensive  set  of  tools  to  simulate  the  performance  of 

OmniGuide  fibers:  hollow,  high  power  transmission  fibers  employing  an  omnidirectional 
dielectric  mirror  as  the  confining  mechanism.  Separately,  we  have  experimentally 
confirmed  the  validity  of  many  of  the  tools.  Our  theoretical/numerieal  analysis  of  field 
propagation  through  OmniGuide  fibers  addressed  several  major  design  aspects  (Fig.  1.1). 
Objectives  for  the  tools  we  developed  were  as  follows: 

1.  Optimal  design  of  OmniGuide  fibers,  including  reduction  of  absorption  and 
radiation  losses. 

2.  Assessment  of  eoupling  efficiency  between  the  high  power  lasers  and  the 
guided  modes  of  an  OmniGuide  fiber. 

3.  Determination  of  beam  quality  degradation  (M^)  due  to  the  inter-modal 
scattering  caused  by  bends  and  other  fiber  imperfections,  including  diameter 
and  ellipticity  variations  along  the  fiber. 

4.  Determination  of  the  thermal  stability  of  the  fibers  including  equilibrium  and 
non-equilibrium  temperature  distributions  across  the  fiber  due  to  CW  and 
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pulsed  sources.  Study  of  the  effects  of  bends  and  other  geometrical  variations 
on  fiber  heating. 

5.  Analysis  of  devices  based  on  OmniGuide  fibers  using  the  Beam  Propagation 
Method 


In  other  efforts,  we  have  analyzed  the  conversion  efficiency  between  the  fundamental 
mode  of  a  low  power  laser  and  the  lowest  loss  TEoi  mode  of  an  OmniGuide  fiber 
(HEn— >TEoi  modal  converters). 


Input 

Coupling  to  HE,  I,  higher  order 
and  cladding  modes 


Region  of 
increased  heating 


Transmission 

Beam  degradation  due  to 
inter-modal  scattering 


Region  of 
increased  heating 


Output 


Beam  quality  degradation 
due  to  higher  order  mode 
content 


Figure  1.1.  Major  areas  of  theoretical  analysis  in  high  power  fiber  design.  Input  - 
coupling  efficiency,  end  face  heating;  Transmission  -  index  profile  design  and 
optimization,  scattering/radiation  due  to  imperfections/bends,  excess  heating  due  to 
bends;  Output  -beam  quality  M^  estimation. 


l.lTooIs  developed  during  Phase  I 

During  phase  I,  the  following  tools  were  developed: 

1)  For  idealized  systems  we  developed: 

a.  Software  to  calculate  the  band  structure  of  a  multilayer  planar  reflector 
using  the  transfer  matrix  method' 

b.  Software  to  calculate  the  band  structure  of  a  cylindrically  symmetrical 
waveguide  with  a  generalized  multilayer  dielectric  structure^’*''®  using 
transfer  matrix  and  Galerkin  type  methods. 
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2)  For  high  index  contrast  systems  slightly  different  (or  perturbed)  from  the 
geometrically  ideal  systems,  we  have  developed  several  computationally  efficient 
perturbation  theories* to  determine  : 

a.  Radiation  absorption  by  the  lossy  dielectric  media. 

b.  Effects  of  waveguide  nonlinearities. 

c.  Propagation  constants  (real  parts  of  the  effective  indexes)  of  the  modes  of 
a  dielectric  waveguide  due  to  an  arbitrary  geometrical  perturbation  of  its 
structure. 

d.  Inter-mode  coupling  due  to  an  arbitrary  geometrical  perturbation  of  a 
waveguide  structure,  including  polarization  mode  dispersion  (PMD). 

e.  Scattering  loss  in  a  hollow  OmniGuide  fiber  due  to  the  roughness  of  the 
reflector’s  surface. 

f.  Changes  in  the  modal  radiation  loss  due  to  the  perturbations  in  reflector’s 
geometry. 

g.  Scattering  losses  due  to  the  roughness  of  the  reflector’s  surface. 

h.  Radiation  along  slightly  non-uniform  waveguides  (propagation-direction 
dependent  perturbation  theory). 

3)  For  OmniGuide  fibers  with  lossy  dielectrics  a  numerical  tool  was  created  based 
on  the  Galerkin  method  to  calculate  complex  propagation  constants.  Furthermore, 
this  software  was  extended  to  calculate  increased  modal  loss  in  bent 
waveguides. 

4)  A  Galerkin-type  numerical  tool  was  developed  to  characterize  coupling  between 
the  modes  of  an  OmniGuide  fiber  and  a  Gaussian  laser  source. 

5)  To  determine  power  capacity,  heat  flow  equations  were  solved  with  a  heat  source 
obtained  from  absorption  of  electromagnetic  energy  during  field  propagation 
through  a  fiber.  Temperature  rise  was  computed  for  CW  and  pulsed  sources,  and 
for  straight  and  bent  fibers. 

6)  For  the  optimization  of  omnidirectional  reflectors,  tools  based  on  nonlinear 
Newton,  conjugate  gradient,  and  non-gradient  based  numerical  optimization 


6 


methods  were  developed  for  the  automated  design  of  planar  multi-reflectors  and 
OmniGuide  fibers  to  obtain: 

a.  Maximum  width  omni-directional  frequency  windows. 

b.  Minimized  absorption  and  radiation  losses. 

c.  Tailored  dispersion  and  reduced  PMD. 

7)  To  analyze  radiation  transmission  along  a  generic  imperfect  waveguide,  software 
implementing  the  Beam  Propagation  Method  (BPM)  was  developed.  Stability  of 
the  software  was  demonstrated  even  for  high  index  contrast  waveguides. 


1.2  Refining  of  theoretical  tools  using  in-house  experimental  data 

Adjustments  and  fine  tuning  of  analytical  models  and  numerical  software  were  made 

throughout  the  program  based  on  extensive  experimental  data  from  our  in-house 

laboratories. 

2  Major  toolboxes  developed  during  phase  I 

Numerical  tools  developed  in  Phase  I  allow  us  to  perform  comprehensive  analyses  of 

all  critical  aspects  of  OmniGuide  fiber  design  and  optimization. 

In  the  following  we  describe  the  major  toolboxes,  their  capabilities,  and  input  and 
output  data  streams.  We  demonstrate  their  utility  by  analyzing  some  typical  OmniGuide 


fiber  designs. 


Toolboxes  develo] 

ped  during  phase  I 

Toolbox  function 

Toolbox  input/output 

Optimal  design  of  OmniGuide  fiber. 
Characterization  of  modal  radiation  and 
absorption  losses.  Optimal  fiber  geometry 
design  for  high  power  applications  with 
reduced  loss. 

Input:  material  parameters  (indexes  of 
refraction  and  bulk  material 
absorption  losses),  core  radius, 
number  of  layers. 

Output:  optimized  individual  layer 
thicknesses  and  corresponding 
losses. 

Laser  to  fiber  coupling  efficiency  analysis. 
Characterization  of  excited  fiber  modes  in  a 
coupler;  mediated  coupling  between  high 
power  laser  and  hollow  Bragg  fibers. 

Input:  coupler  and  waveguide  geometries, 
material  parameters,  diameter  of  the 
incoming  Gaussian  laser  beam. 

Output:  losses  in  the  coupler,  modal 
excitation  coefficients. 

Analysis  of  fiber  performance  degradation 
due  to  geometric  imperfections. 
Characterization  of  increased  modal  loss 

Input:  fiber  geometry,  material  parameters. 
Imperfections  parameters  -  radius  of 
curvature  and  bend  angle  (bends). 
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due  to  generic  planar  bends  and  fiber 
diameter  fluctuations. 

amplitude  of  variations  as  a  function 
of  propagation  length  (diameter 
fluctuations). 

Output:  amplitudes  of  the  excited  modes 
after  the  imperfection  region.  Loss. 

Calculation  of  outgoing  beam  quality  (M  ) 
Outgoing  beam  is  propagated  from  the 
fiber  end  facet  through  free  space  and 
compared  with  Gaussian  beam  divergence. 

Input:  fiber  geometry,  material  parameters, 
amplitudes  of  excited  fiber  modes. 

Output:  electromagnetic  fields  at  the 

specified  distance  from  the  output 
fiber  facet  and  estimates  of  the  M^ 
from  beam  divergence  fit. 

Temperature  distribution  in  hollow  Bragg 
fibers  due  to  CW  and  pulsed  sources. 
Calculation  of  stationary  and  time 
dependent  temperature  rise  along  the  fiber 
for  arbitrary  modal  excitation  pattern. 

Input:  fiber  geometry,  material  parameters, 
amplitudes  of  excited  modes  along 
the  length  of  a  fiber,  CW  or  pulse 
parameters  of  an  incoming  laser. 

Output:  temperature  rise  along  the  fiber  at 
any  point  inside  the  fiber. 

HEii— »TEoi  modal  converter  design. 
Calculates  efficiency  of  conversion  into  the 
lowest  loss  TEoi  mode. 

Input:  fiber  geometry.  Parameters  of  the 
serpentine  mode  converter  -  pitch, 
amplitude,  length  of  the  incoming 
and  outgoing  regions. 

Output:  amplitudes  of  excited  modes. 

Beeim  Propagation  Method.  Modeling  of 
long  distance  radiation  propagation  along 
the  waveguide  in  the  presence  of  slow 
varying  imperfections. 

Input:  fiber  and  slow  imperfections 

geometry  along  the  direction  of 
propagation.  Incoming  profile  of  the 
electromagnetic  fields. 

Output:  profile  of  outgoing  electromagnetic 
fields. 

2.1  Toolbox:  optimal  design  of  OmniGuide  fibers 

Input:  material  parameters  including  indices  of  refraction,  bulk  material  absorption 
losses,  hollow  core  radius,  number  of  mirror  layers. 

Output:  optimized  individual  layer  thicknesses  and  corresponding  fiber  radiation  and 
absorption  losses. 

Using  in-house  transfer  matrix  software  we  first  design  a  planar  multilayer  dielectric 
omnidirectional  reflector  and  then  employ  that  design  to  make  a  cylindrical  internal 
reflecting  surface  to  confine  radiation  in  a  fiber  core.  The  simplest  planar  reflector  design 
employs  a  constant  period  where  thicknesses  of  high  and  low  index  layers  are  fixed 
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across  the  reflector.  Although  the  most  intuitive  and  well-studied,  sueh  reflectors  exhibit 
suboptimal  losses  at  their  wavelength  of  maximum  reflection. 

We  developed  a  design  eode  that  optimizes  the  thieknesses  of  each  individual 
reflector  layer  so  that  the  total  reflection  off  a  planar  stack  is  maximized  for  14  TM,  14  TE 
(HE  II -like)  polarization  of  incoming  radiation.  After  an  initial  optimization  of  the  planar 
reflector,  we  adopt  the  same  layer  thieknesses  for  the  dielectric  mirror  of  an  OmniGuide 
fiber.  Finally,  we  run  another  optimization  code  that  fine-tunes  the  dielectric  layer 
thicknesses  of  the  refleetor  in  the  fiber  geometry. 

The  result  of  the  optimization  process  is  a  dramatic  reduction  in  the  losses  of  the 
OmniGuide  fiber,  compared  to  an  unoptimized  design. 

In  all  the  examples  below  we  will  use  materials**  currently  under  investigation  at 
OmniGuide  Communications  Inc.  (OGCI) 

2.1.1  Example  of  toolbox  application:  optimal  design  of  chalco-polymer  mirror, 
polymer  cladding  system 

In  Progress  Reports  I  and  II  we  considered  a  chalco-polymer-polymer  design  with  the 
following  struetural  parameters.  The  core  radius  Rj  is  chosen  to  be  500  pm.  The  mirror 
materials  are  chosen  to  be  a  chalcogenide  glass  (high  index)  and  a  polymer  (low  index) 
with  indices  of  ni  =  1.5  and  nh  =  2.8,  with  a  bi-layer  thickness,  d=2.S2  pm.  The 
corresponding  material  layer  thicknesses  are  1.97  pm  and  0.84  pm,  respectively.  The 
intrinsic  bulk  absorption  loss  of  the  high  index  material  is  10  dB/m,  an  overestimate  for 
high  purity  chalcogenide  glasses  at  X  =  10.6  pm.  The  cladding  and  low  index  material  are 
assumed  to  have  absorption  loss  of  10*dB/m,  typical  of  polymers.  We  assume  a  mirror  of 
35  layers  (17.5  bi-layers),  thus  Rm  =  549  pm.  At  A,=10.6  pm,  these  structural  parameters 
produee  an  unoptimized fiber  radiation  loss  of  8.9dB/km  (with  a  radiation  loss  decreasing 
by  lOX  with  every  4  bi-layers  added),  and  an  absorption  loss  of  313.6  dB/km. 

Next,  the  first  five  layer  thicknesses  are  optimized  to  reduce  the  absorption  in  the 
polymer  material,  and,  starting  with  the  chalco  layer,  they  are  ehosen  by  the  optimization 
software  to  be  0.51  pm,  0.97  pm,  1.52  pm,  0.97  pm,  and  1.52  pm.  At  X=10.6  pm,  these 
optimized  structural  parameters  produce  a  fiber  radiation  loss  of  3.4  dB/km  (with  a 
radiation  loss  decreasing  by  lOX  with  every  4  bi-layers  added),  and  an  absorption  loss  of 
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79  dB/km.  Thus,  optimization  reduces  the  total  fiber  loss  by  a  factor  of  4  as  demonstrated 
in  Fig.  2.1.1. 


HEjj  dissipation  losses  for  R^=500jjm,  N=35,  nj^/nj=2.8/1.5,  X^p=10.6fj.in 


Figure  2.1.1  Absorption  loss  of  the  chalco-polymer-polymer  fiber  as  a  function  of  the 
wavelength.  For  comparison,  losses  of  the  imoptimized  and  optimized  structures  are 
presented. 

2.2  Toolbox:  coupling  between  high  power  laser  and  OmniGuide  fibers 

Input:  coupler  and  waveguide  geometries,  material  parameters,  diameter  of  the  incoming 

Gaussian  laser  beam.  Coupler  is  modeled  as  a  lossy  metallic  tube  butt-coupled  to 
the  fiber.  Coupler  and  fiber  centers  are  aligned.  Focus  of  the  incoming  Gaussian 
can  be  located  anywhere  along  the  fiber  axis. 

Output:  losses  in  the  coupler,  modal  excitation  coefficients  in  the  fiber. 

In  this  toolbox  we  analyze  coupling  from  a  Gaussian  laser  beam  to  an  OmniGuide 
fiber,  shown  schematically  in  Fig.  2.2.1 .  The  output  from  a  laser  passes  through  a  lens 
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and  is  focused  to  a  spot  of  width  Wq  ,  at  either  the  fiber  endface  or  elsewhere  along  the 

fiber  axis.  From  this  analysis  we  predict  the  amplitudes  of  modes  excited  in  the  fiber. 

The  incident  laser  beam,  assumed  to  have  a  Gaussian  profile,  is  decomposed  into  a 
sum  of  forward-propagating  modes  of  a  homogeneous  waveguide,  and  the  reflected  field 
is  expressed  as  a  sum  of  backward-propagating  modes.  Similarly,  the  field  in  the  fiber  is 
expressed  as  a  sum  over  forward-propagating  fiber  modes.  The  amplitudes  of  these 
modes  are  chosen  so  that  the  electric  and  magnetic  fields  satisfy  the  proper  continuity 
conditions  on  the  fiber  endface.  The  beam  parameters  are  then  varied  to  optimize  the 
coupling  efficiency. 


Figure  2.2.1.  Geometry  of  the  coupling  between  laser  beam  and  an  OmniGuide  fiber.  Due 
to  a  non-  zero  overlap  between  the  incoming  Gaussian  and  the  mirror,  highly  dissipative 
cladding  modes  could  be  excited,  and  an  intermediate  aperture  may  be  necessary. 

2.2.1  Example  of  toolbox  application:  optimal  coupling  to  the  chalco-polymer 
mirror,  polymer  cladding  system 

At  a  wavelength  of  10.6  microns,  we  consider  coupling  to  a  35  layer  OmniGuide  fiber 
with  rifjj  =2.8,  =1.6,  a  core  radius  of  500  microns,  and  a  bilayer  thickness  2.82 

microns.  The  coupling  efficiency  from  a  laser  beam  to  the  HEl  1  core  mode  of  this  fiber 
is  shown  in  Fig.  2.2.2,  below,  as  a  function  of  the  beam  spotsize  .  In  this  example,  the 

beam  was  focused  directly  on  the  fiber  endface  (i.e.,  Zg  =  0).  From  the  figure,  we  see  that 
a  maximum  coupling  efficiency  of  98.6  %  can  be  obtained  from  a  beam  with  Wg=  235 
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Jim.  With  this  beam,  only  0.05  %  of  the  incident  power  is  reflected;  the  remainder  of  the 
incident  beam  excites  undesired  modes  of  the  fiber. 


Spot  Size  (microns) 

Figure  2.2.2.  Coupling  efficiency  from  a  laser  beam  to  the  HEl  1  mode  of  OGCI  fiber  as 
a  function  of  the  spotsize  of  the  beam. 

However,  in  Progress  Report  II  we  established  that  a  more  appropriate  condition  for 
optimal  coupling  is  a  reduction  in  the  excitation  of  the  cladding  (polymer)  modes  (Fig. 
2.2.3),  to  reduce  heating  in  the  polymer  cladding.  This  usually  leads  to  a  suboptimal 
coupling  to  the  HEi  i  mode  at  the  point  of  the  lowest  fiber  end  heating.  Thus,  although  an 
optimal  Gaussian  spot  size  is  235|4.m  for  the  optimal  coupling  to  the  HEn  mode,  coupling 
to  the  cladding  modes  is  still  very  strong  (on  the  order  of  10'^  of  the  incident  power).  In 
the  following  analysis,  we  show  that  power  coupled  to  the  cladding  modes  should  be 
reduced  below  10'^  of  the  incident  power  to  reduce  the  fiber  end  heating  below 
acceptable  values.  For  that,  the  largest  spot  size  we  can  tolerate  without  deliberately 
shielding  the  fiber  end  is  140pm,  at  which  point  the  coupling  to  the  HEn  mode  is  still 
quite  large  and  equal  to  65%.  We  also  established  and  successfully  tested  simple  scaling 
laws  that  approximate  the  coupling  efficiencies  to  the  HEn  and  cladding  modes.  These 
laws  are  derived  assuming  a  Gaussian  HEn  field  profile  of  effective  width  WHEii=235pm, 
and  exponential  decay  of  the  fields  of  an  incoming  Gaussian  laser  beam  in  the  cladding 
region  of  the  fiber.  With  these  assumptions,  we  find  that 
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Figure  2.2.3.  Coupling  to  the  HEl  1  and  cladding  (polymer)  modes  of  the  chalco- 
polymer-polymer  design  as  a  function  of  the  incoming  Gaussian  spot  size. 

2.3  Toolbox:  modeling  increased  modal  loss  due  to  generic  planar  bends  and 
diameter  variations 

Input:  fiber  geometry,  material  parameters,  description  of  imperfection.  The  imperfection 
parameters  for  bending  loss  calculations:  radius  of  curvature  and  bend  angle;  for 
diameter  fluctuations:  amplitude  of  diameter  variation  along  propagation  length. 
Output:  amplitudes  of  the  excited  modes  inside  and  after  the  imperfection  region,  total 
loss  due  to  the  imperfections. 

This  toolbox  models  the  degradation  of  a  field  propagating  through  either  a  bend  or  a 
region  with  diameter  fluctuations.  This  degradation  generally  results  from  scattering  from 
the  incident  fundamental  mode  into  the  higher  order  guided  modes  and  the  radiation 
continuum.  Such  scattering  leads  to  an  additional  loss,  as  well  as  an  increase  in  the  beam 
M^.  We  developed  a  software  tool  based  on  coupled  mode  theory  with  absorbing 
boundary  conditions  that  calculates  fiber  losses  under  such  variations,  and  successfully 
verified  it  against  experimental  data. 
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2.3.1  Example  of  toolbox  application:  bend  losses  in  hollow  metal  fibers  (comparison 
with  experiment)  and  OmniGuide  fibers 

In  Figs.  2.3.1,  2.3.2,  and  2.3.3  we  present  calculations  of  the  straight  and  bend  losses 
of  a  hollow  metallic  (silver)  fiber  with  a  single  dielectric  layer  (silver  iodide),  along  with 
experimental  data  from  Matsuura’®.  In  these  simulations,  the  index  of  refraction  of  silver 
is  taken  to  be  13.5  +  75.3/  (Miyagi  and  Kawakami)  and  the  index  of  refraction  of  the 
silver  iodide  layer  is  2.1.  The  silver  iodide  layer  has  a  thickness  of  0.6  pm.  We  observe 
excellent  agreement  between  the  theoretical  predictions  and  the  experimental  data. 


Figure  2.3.1.  Predicted  and  measured  HEn  mode  loss  of  a  straight  waveguide  as  a 
function  of  fiber  core  radius. 
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Core  Radius  (microns) 

Figure  2.3.2.  Core  radius  =  125  )xm.  Predicted  and  measured  bend  loss  as  a  function  of 
the  bend  curvature.  £  J.  to  plane  of  the  bend.  Measured  for  a  90®  bend. 


Radius  of  Curvature  (cm) 


Figure  2.3.3.  Core  radius  =  160  pm.  Predicted  and  measured  bend  loss  as  a  function  of 
the  bend  curvature.  £  X  to  plane  of  the  bend.  Measured  for  a  90®  bend. 


In  Fig.  2.3.4  we  present  the  bending  losses  after  a  90”  circular  bend  in  the 
unoptimized  OmniGuide  fiber  (described  above)  as  a  function  of  its  core  radius,  for  three 
different  values  of  curvature  radii.  For  the  optimized  design,  bending  losses  improve 
dramatically.  In  Fig.  2.3.5  we  compare  bending  losses  after  a  90”  circular  bend  of  radius 
20cm  between  the  optimized  and  unoptimized  designs. 


Figure  2.3.4.  Core  radius  =160  pm.  Predicted  and  measured  90*^  bend  loss  as  a  function 
of  the  core  radius.  £  ±  to  plane  of  the  bend. 


Figure  2.3.5.  Predicted  total  loss  (linear  plus  bending)  for  chalco-polymer-polymer  OGCI 
fiber  as  a  function  of  the  core  radius.  £  ±  to  plane  of  the  bend.  Simulation  for  a  90” 
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bend,  and  20cm  radius  of  curvature.  The  unoptimised  design  from  Report  I  is  compared 
to  an  optimized  design  presented  in  the  current  report.  Total  loss  is  much  smaller  for  the 
optimized  design,  for  any  value  of  the  core  radius. 


2 

2.4  Toolbox:  calculation  of  spatial  beam  quality  -  M 

Input:  fiber  geometry,  material  parameters,  amplitudes  of  excited  fiber  modes  at  the  fiber 
output  facet. 

Output:  electromagnetic  fields  at  a  specified  distance  from  the  output  fiber  facet  and 
estimates  of  the  M  from  the  beam  divergence. 


Degradation  of  a  field  propagating  through  a  bent  fiber  is  frequently  characterized 
by  the  number.  The  degradation  arises  from  inter-modal  scattering  between  the  initial 
fundamental  mode  and  the  higher  order  modes.  Such  scattering  leads  to  an  increase  in  the 
output  beam  diameter,  an  increase  in  the  angle  of  beam  divergence,  and  “speckling”  of 
the  radiation  profile. 

One  standard  measure  of  the  beam  quality  is  the  number,  which  can  be 
defined  through  the  divergence  of  the  beam  diameter  D(z)  once  light  leaves  the  fiber. 
This  measure  ultimately  reflects  the  smallest  spot  size  that  an  outgoing  beam  can  be 
focused  to.  In  the  Gaussian  beam  approximation,  the  beam  diameter  is  given  by 


Dl,{z)  =  D\0) 


1  + 


?rD\0) 


We  obtain  by  fitting  (z)  to  the  exact  divergence  D(z). 

We  evaluate  the  beam  diameter  D(z)  by  numerically  propagating  the  field  from  the 
output  of  a  fiber  to  a  distance  L  (usually  several  cm)  until  a  reliable  determination  of  M 

is  achieved.  Given  the  field  intensity  I(r,z),  the  beam  diameter  is  defined  as 

1^ 
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The  beam  quality  parameter  is  obtained  by  minimizing  the  discrepancy  Q  between  the 
Gaussian  approximation  and  D(z)  along  the  propagation  length  L; 


L 


D\0) 
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,;rD^(0)j 
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2.4.1  Example  of  toolbox  application:  of  major  propagating  modes 

In  the  following,  we  present  results  of  the  evaluation  of  for  several  low-order 

modes  in  an  OmniGuide  fiber  of  the  chalco-polymer-polymer  design.  The  fiber  is 
composed  of  20  bilayers,  and  has  a  core  ID  of  660pm,  and  an  OD  of  1250pm.  In  Fig. 
2.4.1  we  show  the  calculated  divergence  of  the  TEoi  mode  after  it  exits  the  fiber. 


Figure  2.4.1.  Exact  beam  diameter  D(z)  and  the  beam  optics  approximation  (^)  for 

the  TEoi  mode  of  optimized  fiber.  For  this  mode,  M^=2.17. 

The  following  table  provides  the  numbers  for  some  of  the  main  modes  of  interest 
of  this  OmniGuide  fiber 


HEll 

TEOI 

HE21 

HE12 

EHll 

EH21 

1.10 

2.17 

2.17 

3.02 

3.22 

4.26 
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2.5  Toolbox:  heat  dissipation  in  OmniGuide  fibers  due  to  CW  and  pulsed  sources; 
calculation  of  the  fiber  power  capacity 

Input:  fiber  geometry,  material  parameters,  amplitudes  of  excited  modes  along  the  length 
of  a  fiber,  CW  or  pulsed  source  parameters  of  the  incoming  laser  beam. 

Output:  temperature  rise  along  the  fiber  at  any  point  inside  the  fiber.  Possible  fiber 
configurations  include  straight,  straight  with  a  coupler,  bent. 

This  toolbox  calculates  the  stationary  temperature  distribution  along  the  fiber  due  to  a 
CW  source,  and  the  peak  temperature  distribution  along  the  fiber  due  to  a  pulsed  laser 
source.  We  developed  a  rigorous  solver  of  the  heat-transfer  equations  with  exact  heat 
sources  due  to  absorption  of  electromagnetic  energy  in  the  fiber  material.  This  solver  also 
allows  investigation  of  heating  due  to  imperfect  laser-fiber  coupling,  and  fiber  bends. 

In  general,  an  OmniGuide  fiber  consists  of  a  hollow  core  surrounded  by  alternating 
layers  of  materials  with  high  and  low  indices  of  refraction.  Fig.  2.5.1  (a)  illustrates  the 
cross  section  of  the  fiber.  Light  is  guided  primarily  in  the  hollow  core  of  radius  i?,-  by 
means  of  the  dielectric  mirror  having  N  dielectric  bi-layers  each  of  thickness  d  extending 
to  the  radius  R^-  Outside  of  the  mirror  there  is  a  cladding  of  radius  Rg. 

The  mirror  is  designed  in  such  a  way  that,  for  a  given  frequency,  the  modal  fields 
decay  exponentially  in  the  mirror  region.  Thus,  field  penetration  into  the  mirror  is  limited 
largely  to  the  first  few  bi-layers.  Moreover,  as  the  core  area  is  100-1000  times  larger  than 
the  mirror  area,  power  dissipation  in  the  mirror  due  to  the  material  absorption  is  greatly 
suppressed.  In  the  following,  we  refer  to  the  power  dissipated  in  the  first  few  bi-layers  of 
the  mirror  due  to  the  finite  modal  field  penetration  as  Hats  (Fig.  2.5.1  b)). 

Another  source  of  power  dissipation  is  absorption  of  the  radiation  leaking  from  the 
finite-sized  mirror.  This  radiation  will  be  mostly  converted  to  heat  in  the  over  clad  region. 
The  radiated  power  dissipated  in  this  region  is  defined  as  Hrad-  Fibers  can  be  built  that  can 
withstand  kilowatts  of  transmitted  power,  even  without  forced  cooling,  by  designing 
them  so  that  most  of  the  energy  propagates  in  the  hollow  core,  minimal  energy 
propagates  in  the  mirror,  and  minimal  radiation  leaks  into  the  over  clad. 
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Figure  2.5.1.  Power  capacity  analysis  of  a  hollow  Bragg  fiber,  (a)  Schematic  of  a  fiber. 
Hollow  core  of  radius  Rj,  mirror  region  of  N  bilayers  each  of  thickness  d  surrounded  by 
an  over-cladding  extending  to  Rq.  (b)  Two  major  types  of  heat  sources  (shown  in  light 
gray)  -  Habs  is  due  to  a  field  penetration  and  absorption  in  the  mirror  region,  Hrad  is  due  to 
the  radiation  leakage  through  a  finite  size  mirror  and  absorption  in  the  over-cladding. 


2.5.1  Example  of  toolbox  application:  fiber  heating  with  CW  laser  source  (general) 

We  first  present  an  analysis  of  the  temperature  increase  in  a  fiber  transmitting  P  [W] 

of  power  coming  from  a  CW  laser  source.  Cooling  is  provided  by  either  dry  air  or  water 
convection  flow  at  room  temperature,  Tc  [°C].  Heat  transfer  fi'om  fiber  to  coolant  is 
characterized  by  a  surface  heat-transfer  coefficient  h  [W/m^/K].  For  CW  sources,  the 
temperature  rise  in  a  fiber  is  very  sensitive  to  the  cooling  conditions,  with  most  of  the 
temperature  rise  occurring  at  the  cladding-air/water  interface,  and  only  a  small 
temperature  difference  is  observed  across  the  fiber  profile. 

We  assume  laser  beam  coupling  to  a  set  of  fiber  modes  with  radial  index  i  and 
angular  index  m,  with  amplitudes  C,_,„(2)  along  the  length  of  a  fiber.  Due  to  modal 

interference,  the  hottest  region  in  the  fiber  is  not  necessarily  at  the  coupling  facet,  but  is 
determined  by  a  particular  modal  excitation  pattern  and  is  usually  somewhere  near  the 
coupling  facet.  The  equilibrium  temperature  in  this  region  can  be  found  by  solving  the 
time-independent  heat  transfer  equation  with  insulating  boundary  conditions  on  the  fiber 


end  facets.  The  heat  source  is  computed  exactly  from  the  excited  modal  pattern.  Energy 
loss  in  the  layers  with  high  absorption  will  contribute  most  to  the  heating;  therefore  the 
heat  source  will  be  a  very  non-uniform  function  in  the  radial  direction  (Fig.  2.5.2). 

Figure  2.5.2.  Modeling  of  fiber  heating.  Most  of  heating  comes  from  energy  dissipation 
in  the  mirror  layers  made  of  highly  absorbing  material  (in  red).  This  dissipation  defines  a 

Jo=h(T(Ro)-Tamb) 


k  -  mirror  heat  conductance 


heat  source  very  non-uniform  in  the  radial  direction.  Heat  conductance  to  the  core  and 
through  the  fiber  end  facets  is  assumed  zero,  with  all  the  heat  assumed  to  dissipate 
through  the  outer  surface  of  the  fiber. 

In  Fig.  2.5.3  we  demonstrate,  for  example,  the  angle-averaged  heat  source  arising 
from  dissipation  of  the  HEn  mode  of  a  large  core  (Ri=330pm)  OmniGuide  chalco- 
polymer-polymer  fiber  at  z  =  0.  As  we  pointed  out  before,  most  of  heating  comes  from 
energy  dissipation  in  the  mirror  layers  made  of  the  more  highly  absorbing  material  -  the 
polymer.  This  dissipation  creates  a  heat  source  that  is  very  non-uniform  in  radius,  with 
maxima  in  the  highly  absorbing  layers. 
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Figure  2.5.3.  Heat  source  S„{p,z)\^^^  of  the  HEn  mode  in  the  fiber  mirror  region.  Radius 

of  the  fiber  core  is  Ri=330|j,m.  All  the  heating  comes  from  energy  dissipation  in  the 
mirror  layers  made  of  highly  absorbing  material.  This  dissipation  defines  a  very  non- 
uniform  heat  source  in  radius  with  maxima  in  the  highly  absorbing  layers.  Exponential 
decay  of  the  modal  field  inside  of  the  reflector  is  responsible  for  most  of  the  heat 
dissipated  in  the  first  few  mirror  layers. 

In  Progress  Reports  II  and  III  we  presented  a  rigorous  solution  of  the  time- 
independent  heat  flow  equation  with  a  heat  source  arising  from  absorption  of  the 
propagating  modes.  We  also  established  that  if  interference  effects  can  be  neglected  the 
peak  temperature  value  Ti  is  occurs  at  the  inner  core  of  the  fiber  end  facet,  and  is  given 
by 
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2.5.2  Example  of  toolbox  application:  increased  fiber  heating  with  CW  laser  source 
due  to  imperfect  coupling  to  the  HEn  mode 

In  the  following,  we  analyze  the  heating  of  the  fiber  when  laser  light  is  coupled 

through  a  coupler  that  excites  higher  order  modes.  The  OmniGuide  chalco-polymer- 
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polymer  fiber  in  this  example  was  designed  for  optimal  performance  at  10.6pm  with  a 
fiber  core  radius  of  Ri=330pm.  A  typical  tube  coupler  (Fig.  2.5.4)  used  in  our 
experiments  consists  of  a  short  metallic  tube  inserted  into  the  fiber’s  hollow  core  through 
a  cooled  metallic  block. 


Figure  2.5.4.  High  power  coupling  to  the  fiber  via  tube  coupler.  Tube  coupler  (in  blue)  is 
inserted  into  the  block  of  metal  (also  in  blue)  to  shield  fiber  end  from  direct  illumination 
by  laser  beam.  Coupler  tube  is  a  thin  metallic  tube  of  inner  core  radius  close  to  that  of  the 
fiber  core  radius. 

Depending  upon  the  coupling  conditions,  such  as  the  diameter  of  the  incoming 
Gaussian  laser  beam,  higher  order  modes  in  the  fiber  can  be  excited.  For  the  particular 
fiber/coupler  geometry  used  in  our  experiments,  we  calculated  that  the  incident  power 
was  distributed  in  the  following  way;  80%  in  the  HEn  mode,  10%  in  the  HEn  mode  and 
the  rest  in  higher-order  modes.  The  losses  of  the  HEn  and  HE12  modes  are  0.2  dB/m  and 
1  dB/m,  respectively.  Although  the  excitation  of  the  HEn  mode  is  eight  times  larger  than 
the  HE  12  mode,  its  losses  are  five  times  smaller  then  the  losses  of  the  higher  order  mode. 
Therefore  the  heat  sources  produced  by  dissipation  of  these  two  modes  are  comparable  in 
magnitude.  Moreover,  because  of  constructive  and  destructive  interference  between  the 
modes,  the  losses  can  vary  substeintially  along  the  length  of  the  fiber,  leading  to  the 
appearance  of  “hot”  spots.  In  Fig.  2.5.5  we  plot  temperature  rise  along  3m  of  fiber, 
assuming  an  incident  power  of  1 W,  and  a  surface  heat  transfer  coefficient  of  h  =  5 
[W/m^/K].  Note  that  the  HE]  1  mode  alone  contributes  2.4  /W  rise  in  temperature 
almost  uniformly  across  the  length  of  the  fiber  (green  curve  on  Fig.  2.5.5),  while  other 
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excited  modes  have  a  larger  contribution  to  the  heating.  Also,  the  temperature  distribution 
along  the  fiber  length  is  very  nonuniform  (red  curve),  with  hot  spots  appearing  about 
every  3cm,  corresponding  to  the  beat  length  between  the  HEn  and  HE12  modes.  Note  also 
that  if  interference  effects  between  modes  are  not  taken  into  account,  the  heating  profile 
(blue  curve)  only  describes  the  average  heat  rise  along  the  fiber  length,  totally  missing 
the  phenomenon  of  “hot”  spots  -  an  important  failure  mechanism  in  high-power  fibers. 


Rc^^=33Ci|im,  R^^^j^=625nm.  X=10.6|im 


Figure  2.5.5.  Temperature  rise  per  1 W  of  input  power  along  a  3m  fiber  with  heat  transfer 
coefficient  h=5  [W/m/K].  Green  curve  shows  temperature  rise  when  only  HEn  mode  is 
excited.  Red  curve  shows  temperature  rise  when  all  the  modes  are  excited  with  80%  in 
HEn  mode  and  10%  in  HE12  mode.  Strongly  oscillatory  nature  of  the  temperature 
distribution  is  due  to  interference  between  excited  fiber  modes.  Blue  curve  corresponds  to 
the  temperature  distribution  when  modal  interference  is  not  taken  into  accoxmt. 

These  predictions  were  tested  experimentally  by  taking  a  thermal  image  (Fig. 
2.5.6)  of  a  short  fiber  piece  near  the  tube  coupler.  As  predicted,  we  see  “hot”  spots 
appearing  every  2-5cm  along  the  fiber  length.  A  more  rigorous  comparison  between 
theory  and  experiment  is  ongoing. 
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Figure  2.5.6.  Thermal  image  of  the  50cm  fiber  piece  near  the  tube  coupler  with  several 
Watts  of  incoming  power.  Note  the  temperature  oscillations  every  2-5cm  appearing  as 
“hot”  spots  along  the  fiber  length.  Bright  light  spot  in  the  upper  left  comer  is  a  reflection 
from  the  lens.  Cube-shaped  object  close  to  the  lens  is  a  coupler.  Bright  circular  objects 
are  cooling  fans  with  their  hot  motors  glowing  in  the  IR. 


2.5.3  Example  of  toolbox  application:  increased  fiiber  heating  with  CW  laser  source 
due  to  a  fiber  bend 

When  bends  are  present  in  the  fiber,  mode  conversion  is  induced.  The  amount  of 
this  mode  conversion  increases  with  the  fiber  core  radius.  Higher  order  modes  excited  by 
the  bend  have  much  larger  attenuation,  leading  to  an  increase  in  heating  in  the  fiber  bend. 
In  Fig.  2.5.7,  we  show  the  temperature  rise  along  a  90°  circular  bend  with  a  20cm  bend 
radius.  We  assume  an  HEn  mode  carrying  1  Watt  of  power  is  incident  on  the  bend.  At 
the  beginning  of  the  bend,  before  substantial  mode  conversion  has  occurred,  heating  is 
dominated  by  the  HEi  i  modal  loss,  leading  to  a  temperature  rise  of  2.4  °C/W.  Further  into 
the  bend,  modal  conversion  is  substantial  leading  to  an  average  temperature  rise  of  1 8 
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°C/W.  Note  the  appearance  of  the  “hot”  spots  along  the  bend  length  because  of  the  modal 
interference. 


90"bendofR=20cm,R  =330^111,  R  ,  =625nni,N,  =18,  X=10.6|im 

’  core  ^  ’  cooler  ^  ^  layers  ^ 


Figure  2.5.7.  Heating  along  90°  fiber  bend  at  20cm  bend  radius,  assuming  excitation  by 
HEii  mode  with  1  Watt  of  power.  The  heat  transfer  coefficient  h=5  [W/m/K].  At  the 
begirming  of  the  bend,  before  substantial  mode  conversion,  heating  is  dominated  by  the 
HEi  I  modal  loss  leading  to  a  temperature  rise  of  2.4  °C/W.  After  2cm  into  the  bend 
modal  conversion  is  substantial  leading  to  an  average  temperature  rise  of  1 8  °CfW. 


2.5.4  Example  of  toolbox  application:  increased  fiber  heating  with  a  pulsed  laser 
source  (general) 

We  now  analyze  the  thermal  behavior  of  an  OmniGuide  fiber  with  a  pulsed  power 
source.  We  assume  Ppeak  [W]  peak  power  is  transmitted  through  the  fiber  with  a  pulse 
duration  t  [s]  and  a  repetition  rate  T  [s].  We  define  Pav  to  be  the  average  transmitted 

power:  Pav=PpeakT/T.  Commonly,  the  pulse  duration  is  much  smaller  than  the 

2 

characteristic  temperature  equilibration  time  across  the  fiber  cross-section  Xeq^CvpR  /k. 
In  this  short  pulse  regime,  the  temperature  distribution  across  the  fiber  cross-section  will 
be  modified  from  that  of  the  steady-state  distribution  corresponding  to  the  average 
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transmitted  power  Pay.  In  general,  power  handling  of  a  pulsed  source  is  more  challenging 
than  power  handling  of  CW  radiation  of  the  same  average  intensity.  Physically,  this  can 
be  attributed  to  the  excitation  by  a  time-dependent  source  of  localized  “heat  waves” 
(Figure  2.5.8). 


Figure  2.5.8.  Phenomenon  of  “heat  wave”.  Heat  dissipation  from  a  temporally  periodic 
heat  source  deteriorates  as  temporal  frequency  of  the  source  gets  larger,  leading  to  an 
increase  in  local  heating.  Heat  from  a  time  dependent  source  is  localized  inside  of  a 
region  of  size  L,  inversely  proportional  to  the  square  root  of  the  source  frequency. 


The  larger  the  frequency  of  the  heat  somce,  (i.e.,  the  shorter  the  pulse  duration)  the 
more  localized  is  the  “heat  wave”,  thus  weakening  heat  exchange  through  the  outer 
boundary  of  the  fiber  to  the  surrounding  coolant.  For  the  very  lossy  cladding  modes, 
heating  becomes  even  more  problematic  because  of  the  additional  localization  of  the  heat 
in  the  direction  of  propagation  due  to  the  very  rapid  modal  decay.  Special  care  should  be 
taken  to  reduce  coupling  to  the  cladding  modes. 

As  in  the  analysis  of  CW  sources,  we  assume  laser  beam  coupling  to  a  set  of  fiber 
modes  with  radial  index  /,  angular  index  m,  and  frequency  ©,  with  excitation  coefficients 
C,  „,(2)  along  the  length  of  the  fiber.  The  steady-state  temperature  in  the  fiber  can  be 

found  by  solving  the  time-dependent  heat  transfer  equation  with  boundary  condition  of 
zero  flux  across  the  fiber  side  facets.  The  heat  source  is  computed  exactly  from  the 
excited  modal  pattern. 
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In  Progress  Report  III,  we  presented  a  rigorous  solution  of  the  time-dependent  heat 
flow  equations  with  a  heat  source  arising  from  absorption  of  the  propagating  modes. 


2.5.5  Example  of  toolbox  application:  heating  due  to  pulsed  sources  in  chalco- 
polymer-polymer  OGCI  fibers 

In  Progress  Reports  II  and  III,  we  established  that,  neglecting  modal  interference, 
the  hottest  region  of  the  fiber  is  on  the  coupling  facet  of  the  fiber  at  the  core/mirror 
boundary.  The  peak  temperature  value  Tj  in  this  region  over  time  was  approximated  by 
solving  a  simplified  time-dependent  heat  transfer  equation  with  boundary  conditions  of 
zero  flux  across  the  coupler/fiber  facet  and  the  end  facet,  and  by  assuming  that  all  the 
heat  flux  was  coming  into  the  fiber  from  a  very  narrow  heat  source  at  the  fiber 
core/mirror  boundary.  The  maximum  temperature  was  then  established  to  be 
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Thus,  assuming  that  high-loss  cladding  modes  » 1  are  suppressed,  the 


temperature  rise  for  the  pulsed  source  as  a  function  of  the  average  power,  pulse  duration 

1  T -T 

and  pulse  repetition  rate  will  be  given  by  a  universal  curve  in  coordinates  ( — -j=,— - 

Pa, 

In  these  coordinates,  the  temperature  rise  is  given  by  a  family  of  straight  lines  each 
corresponding  to  a  different  cooling  mechanism.  The  CW  power  capacity  is  calculated  by 

letting — ^  =  0 .  For  example,  a  chalco-polymer  design  with  RpSSOpm  and  14  bilayers, 
Vyjr 

a  maximum  allowed  temperature  rise  of  AT=150°C,  dry  air  cooling,  with  HEn  mode 
propagation,  has  a  CW  power  capacity  of  Pav=73W  .  The  CW  temperature  rise  of  this 
fiber  will  be  Tcw=2”C/W.  For  a  pulsed  cutting  laser  with  T=200ns  and  v=100Hz,  the 
maximum  power  capacity  is  reduced  to  Pav=62W,  with  Ppeak=3  lO^W. 


28 


The  rigorous  calculation  of  thermal  behavior  with  pulsed  sources  described  in 
Progress  Report  III  shows  that  the  temperature  rise  is  actually  smaller  when  evaluated 
with  a  realistic  heat  source,  rather  than  a  6-function-like  heat  source.  For  example,  for  the 
same  chalco-polymer-polymer  design,  as  the  repetition  rate  of  the  pulses  increases,  the 
temperature  distribution  tends  to  get  closer  to  the  CW  distribution  with  the  same  average 
power.  In  Fig.2.5.9  we  present  the  temperature  rise  due  to  a  pulsed  source  as  a  function  of 

the  pulse  parameter  — The  simplified  theory  with  a  6-function-like  heat  source 

predicts  a  straight  line.  Detailed  calculations  show  that  this  line  gives  an  upper  bound  on 
the  actual  heating  due  to  a  pulsed  source.  In  reality,  because  of  the  finite  extent  of  the 
heat  source,  the  temperature  rise  curves  tend  to  a  CW  temperature  rise  as  the  pulse 
repetition  rate  increases.  In  Fig.2.5.9,  several  temperature  rise  curves  are  presented  as  a 

function  of  .  They  are  labeled  and  calculated  for  particular  pulse  repetition  rates  v, 
where  the  pulse  duration  x  is  varied. 


Figure  2.5.9.  Temperature  rise  curves  as  a  function  of  the  laser  pulse  parameter;  v  is  the 
pulse  repetition  rate,  x  is  the  pulse  duration.  Different  dashed  curves  are  calculated  for  a 
realistic  heat  source,  and  correspond  to  different  pulse  repetition  rates.  Solid  blue  curve  is 
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calculated  for  a  simplistic  6-flmction  like  heat  source  which  provides  an  upper  bound  on 
the  heating  increase  due  to  pulsed  sources. 

2.6  Toolbox:  HEn— >TEoi  modal  conversion 

Input:  fiber  geometry,  parameters  of  a  serpentine  mode  converter.  Serpentine  mode 

converter  is  characterized  by  the  pitch,  amplitude,  and  length  of  the  incoming  and 
outgoing  transition  regions. 

Output:  amplitudes  of  excited  modes. 

As  described  in  detail  in  Progress  Report  I,  the  HEn  and  TEoi  modes  are  the  most 
desirable  for  operation  in  an  OnmiGuide  fiber.  The  HEn  mode  has  the  best  coupling  to  a 
Gaussian  laser  beam,  while  the  TEoi  mode  has  the  lowest  transmission  loss,  but, 
unfortunately,  does  not  couple  well  to  an  incident  laser  beam.  With  a  low-loss 
HEii-^  TEoi  mode  converter,  however,  the  HE] i  mode  can  be  employed  for  efficient 
coupling,  and  the  TEoi  mode  can  be  used  for  optimal  transmission.  We  have  successfully 
initiated  a  theoretical  and  experimental  effort  to  develop  a  mode  converter  that  efficiently 
performs  this  conversion  in  an  OnmiGuide  fiber. 


2.6.1  Example  of  toolbox  application:  HEn  (LPoi)— ^TEoi  (LPn)  modal  conversion  in 
SMF-28  fiber  at  3i=633nm 

A  prototype  mode  converter  that  we  studied  consists  of  a  deformation  of  the  fiber  axis 
into  an  appropriate  serpentine  bend.  This  deformation  is  achieved  by  pressing  the  fiber 
between  two  sets  of  wires  wrapped  on  two  cylinders  (Figure  2.6.1).  The  wire  diameter  is 
chosen  so  that  the  bend  spacing  in  the  fiber  axis  matches  the  beat  length  between  the 
HEn  and  TEoi  modes.  Initially,  all  the  power  is  launched  into  the  HEn  mode.  If  the 
diameter  of  the  wire  matches  exactly  the  beat  length  between  the  modes,  almost  100% 
conversion  efficiency  can  be  achieved  by  varying  the  pressure  on  the  top  set  of  wires  to 
induce  the  proper  deformation.  In  practice,  this  matching  condition  is  hard  to  achieve,  so 
conversion  efficiency  is  usually  smaller. 
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Figure  2.6.1.  Mode  converter  induces  a  serpentine  2D  bend  by  pressing  the  fiber 
between  two  sets  of  wires  wrapped  on  two  cylinders.  The  diameter  of  wires  is  matched  to 
the  beat  length  between  two  modes.  Amplitude  of  fiber  oscillations  6  is  chosen  to 
optimize  mode  conversion. 

In  Figure  2.6.2  we  present  simulations  of  the  maximum  converter  efficiency  based  on 
a  specialty  step-index  silica  fiber  operating  at  X,=633nm.  At  this  wavelength,  the  specialty 
fiber  has  a  single  HEn  (LPoi)  mode  and  three  almost  degenerate  LPn  modes,  one  of 
which  is  the  TEoi,  above  cutoff  The  fiber  has  a  core  radius  of  Rc=2.23  pm,  and  an  index 
contrast  An/nc=0.34%.  The  beat  length  between  the  HEn  and  TEoi  modes  for  this  fiber  is 
247  pm.  In  the  experimental  realization  of  the  mode  converter,  14  consecutive  wire 
wraps  of  diameter  D=254  pm  were  employed. 


Rxod  ends.  6«0.lM.m,V633nm,  specially  flben  d2«i|im 


Figure  2.6.2.  Mode  conversion  in  a  serpentine  2D  bend  along  the  length  of  a  converter. 
Polarization  of  the  initial  HEi  i  mode  corresponds  to  £  i.  to  plane  of  the  bend. 
Maximum  58%  mode  conversion  into  LPn  is  predicted  with  29%  ending  up  in  the  TEoi 
mode  for  the  optimal  amplitude  of  fiber  bend  oscillation  of  6=100nm. 
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We  used  in-house  software  to  extract  the  modal  content  from  intensity  images  of  the 
converter  output  that  were  recorded  by  a  CCD  camera.  From  the  experimental  images, 
(Figure  2.6.3)  we  found  that  we  achieved  55%  conversion  into  the  LPn  group,  which 
agreed  very  well  with  the  predicted  conversion  rate. 


Figure  2.6.3.  CCD  images  of  mode  conversion  in  a  serpentine  2D  bend.  Left  -  HEn 
mode  in  initial  unperturbed  SMF-28  fiber.  Right  -  mode  conversion  by  gently  pressing 
on  the  top  set  of  wires  -  maximum  conversion  of  58%  into  LPu  achieved  by  varying  the 
pressure  on  the  top  set  of  wires. 


2.7  Toolbox:  Beam  Propagation  Method 

Input:  description  of  the  fiber,  description  of  (gentle)  imperfections  along  the  direction  of 
propagation.  The  incoming  profile  of  the  electromagnetic  field  has  to  be  provided 
to  start  the  propagation.  The  BPM  software  has  its  own  interface  based  on  the 
programming  language  Guile. 

Output:  electromagnetic  fields  along  the  direction  of  propagation. 

The  general  problem  of  field  propagation  through  an  index  profile  with  arbitrary 
(but  slowly-varying)  perturbations  can  be  in  many  cases  be  addressed  by  the  Beam 
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Propagation  Method  (BPM).  Currently,  all  finite-difference  BPM  codes  of  which  we  are 
aware,  including  commercially  available  codes,  are  unstable  when  applied  to  the  study  of 
high  index  contrast  systems  -  the  computed  fields  grow  arbitrarily  large  during 
propagation.  While  there  is  some  evidence  of  stable  Finite  Element  BPM  codes,  they  are 
not  yet  in  the  public  domain  and  are  definitely  not  optimized  for  the  study  of  high  index 
contrast,  multi-layered  index  profiles.  We  have  developed  BPM  software  that  is 
numerically  stable  and  includes  a  wide  angle  correction  based  on  and  2"'*  order  Pade 
expansions.  This  toolbox  allows  us  to  study  complex  devices  built  out  of  OmniGuide 
fibers,  as  well  as  analyze  the  effects  of  imperfections  present  in  these  fibers. 


2.7.1  Example  of  toolbox  application:  field  propagation  in  fiber  with  microbending 

A  common  imperfection  in  a  fiber  is  microbending.  We  present  simulation  results 

of  TEoi  mode  propagation  through  an  S-shaped  microbend  with  bending  radius 
Rbend-500p.m  and  bending  angle  of  12°.  The  OmniGuide  fiber  has  a  core  radius  of  2.5|im 
and  7  mirror  layers  optimized  for  transmission  at  1.55p,m.  In  Fig.  2.7.1  we  show 
snapshots  of  the  TEoi  energy  density  in  the  fiber  cross-section  along  the  bend.  The  BPM 
calculation  predicts  a  2dB  loss  from  the  bend. 


Figure  2.7.1.  Energy  dissipation  of  TEoi  mode  in  an  S-bend  with  Rbend=I  000pm  and 
bending  angle  of  12°.  A  sequence  of  energy  density  snapshots  in  the  fiber  cross  section 
along  the  fiber  is  presented.  Outmost  left  picture  corresponds  to  the  starting  energy 
density  before  the  bend.  Outmost  right  picture  corresponds  to  energy  density  after  the 
bend.  Colorscale  is  red-maxumum  energy  density,  white-minimum. 
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2.7.2  Stability  analysis  of  BPM  code  (The  work  presented  in  this  section  was 
performed  by  Stanford  University) 

The  effect  of  the  Crank-Nicolson  a  parameter  on  the  numerical  stability  and 
dissipation  of  the  OGBPM  simulation  code  has  been  investigated.  Field  propagation  in 
the  OGBPM  is  accomplished  with  the  Crank-Nicolson  method,  and  it  is  employed  as 
follows: 

^  =  d(z)f(2) 

02 

F{z  +  dz)  =  (1  -  dzad(z  +  dz))(}  +  dz(\  -  a)d{z))F{z) 
where  dz  is  the  step  size  in  the  direction  of  propagation  through  which  the  field  F  is 
advanced. 

In  the  following,  we  investigate  the  performance  of  the  OGBPM  code  as  a  function  of 
the  Crank-Nicolson  parameter  a.  Previously,  we  observed  that  stability  of  OGBPM  can 
be  substantially  improved  by  increasing  a  from  its  ideal  value  of  V2,  where  the  Crank- 
Nicolson  algorithm  is  second-order  accurate.  However,  when  a  is  changed  from  its  ideal 
value,  numerical  dissipation  is  introduced  into  the  system.  Moreover,  numerical 
convergence  also  suffers,  with  the  error  of  the  method  now  being  proportional  to  dz 
rather  than  dz^  for  a  =  1/2.  The  purpose  of  this  report  is  to  investigate  these  tradeoffs 
between  numerical  stability  and  munerically  induced  dissipation. 

To  investigate  the  effect  of  the  a  parameter,  an  index  guided  fiber  is  considered 
(Fig.4.1). 


Figure  2.7.2  4a  x4a  computational  cell,  dielectric  cylinder  (£  =  1 1.56)  with  diameter  a,  in 
vacuum. 
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The  initial  field  profile  is  imported  from  the  MIT  Photonic  Band  (MPB)  simulation 
code  with  the  same  fiber  structure.  The  frequency  is  set  to  the  eigen-frequency  of  the 
mode  (0.357954).  The  wavevector  is  set  to  1  as  in  the  MPB  simulation.  The  following  are 
typical  distance  vs.  total  energy  plots  after  propagation  of  the  initial  field  by  OGBPM 
(Fig.  2.7.3) 
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Figure  2.7.3.  Distance  vs.  total  energy  plots  after  propagation  of  initial  field  by  OGBPM. 
Diverging  (a  ==  0.55)  and  non-diverging  (a  =  0.70)  examples. 


The  initial  drop  in  energy  is  possibly  due  to  a  mode  mismatch  between  the  MPB  code 
and  OGBPM.  One  can  fit  the  decay  (growth)  in  both  cases  with  exponential  functions  of 
the  form  A*exp(Bz).  In  both  cases,  a  very  good  fit  is  achievable,  indicating  that  the  decay 
(growth)  is  typical  numerical  divergence  /  dissipation  (Fig.  2.7.4). 
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Figure  2.7.4.  Fit  of  the  energy  decay  curves  with  exponential  functions  of  the  form 
A*exp(Bz). 


The  relationship  between  a  and  the  growth/decay  constant  B  is  summarized  in  the 
tables  below.  Observe  that  a=0.62  provides  stable  propagation  with  no  dissipation. 
Choosing  (x>0.62  leads  to  a  stable  propagation  scheme,  but  an  artificial  numerical 


dissipation  is  introduced.  Choosing  a<0.62  leads  to  an  unstable  propagation  scheme. 
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chalcogenide  glass  for  the  low  index  layer  and  the  other  uses  a  polymer. 
Properties  of  the  materials  are  as  follows: 


Polymer 

Chaleo  Glass 

Wavelength 

X  [fim] 

2-11 

2-11 

Index  of  refraction 

N 

1.5 

2.3-2.8 

Bulk  material  absorption 

a  [dB/m] 

10’ 

1  0 

Heat  conductivity 

A:[W/m/K] 

0.16-0.27 

0.252 

Density 

p[kg/m'] 

1400 

4400 

Specific  heat 

Cv[JAg4C] 

1000 

294 
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